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Abstract 

A model with six mutually-exclusive compartments related to Dengue disease is 
presented. In this model there are three vector control tools: insecticides (larvicide and 
adulticide) and mechanical control. The problem is studied using an Optimal Control 
(OC) approach. The human data for the model is based on the Cape Verde Dengue 
outbreak. Some control measures are simulated and their consequences analyzed. 

Key words: modeling, optimal control, basic reproduction number, vector control 
MSG 2000: 92B05, 93C95. 

1 Introduction 

Dengue is a vector borne disease transmitted to humans by the bite of an infected female 
Aedes mosquito. Dengue transcends international borders and can be found in tropical 
and sub-tropical regions around the world, predominantly in urban and semi-urban areas. 
The risk may be aggravated further due to climate changes and to the globalization, as a 
consequence of the huge volume of international tourism and trade [10J. 

There are four distinct, but closely related, viruses that cause Dengue. Recovery from 
infection by one virus provides lifelong immunity against that virus but confers only par- 
tial and transient protection against subsequent infection by the other three viruses 
Unfortunately, there is no specific effective treatment for Dengue. 
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Primary prevention of Dengue resides mainly in mosquito control. There are two pri- 
mary methods: larval control and adult mosquito control, depending on the intended target 
[6]. The application of adulticides can have a powerful impact on the abundance of adult 
mosquito vector. However, the efficacy is often constrained by the difficulty in achieving 
sufficiently high coverage of resting surfaces. This is the most common measure. Besides, 
the long term use of adulticide has several risks: the resistance of the mosquito to the 
product reducing its efficacy and the killing of other species that live in the same habitat. 
Larvicide treatment is done through long-lasting chemical in order to kill larvae and prefer- 
ably have WHO clearance for use in drinking water [1]. Larvicide treatment is an effective 
way to control the vector larvae, together with mechanical control, which is related with 
educational campaigns. The mechanical control must be done both by public health officials 
and by residents in affected areas. The participation of the entire population is essential to 
remove still water from domestic recipients, eliminating possible breeding sites |12j . 

Mathematical modeling is an interesting tool for understanding epidemiological diseases 
and for proposing effective strategies to fight them [5]. 

2 Mathematical Model 

Taking into account the model presented in [3] and the considerations of [8], it is 
proposed a new model more adapted to the Dengue reality. The notation used in the 
mathematical model includes three epidemiological states for humans, indexed by h: 

Sh susceptible (individuals who can contract the disease); 

Ih infected (individuals capable of transmitting the disease to others) ; 

Rh resistent (individuals who have acquired immunity). 

It is assumed that the total human population (Nh) is constant and Nh = Sh(t) + 
Ih{t) + Rh(t) at an Y time t. 

There are three other state variables, related to the female mosquitoes, indexed by m: 

A m aquatic phase (that includes the egg, larva and pupa stages); 

S m susceptible (mosquitoes that are able to contract the disease); 

I m infected (mosquitoes capable of transmitting the disease to humans). 

Due to short lifetime of mosquitoes (approximately 10 days), there is no resistant phase. 
Humans and mosquitoes are assumed to be born susceptible. 

To analyze the effect of campaigns in the combat of the mosquito, three controls are 
considered: 

ca proportion of larvicide, and < ca < 1; 

c m proportion of adulticide, and < c m < 1; 

a proportion of mechanical control, and < a < 1. 
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The aim of this work is to simulate different realities in order to find the best policy 
to decrease the number of human infected. A temporal mathematical model is introduced, 
with mutually-exclusive compartments, to study the outbreak of 2009 in Cape Verde islands 
and improving the model described in [TJ. 

The model considers the following parameters: 

Nh total population; 

B average daily biting (per day); 

/3 m h transmission probability from I m (per bite); 

/3hm transmission probability from Ih (per bite); 

1/ Hh average lifespan of humans (in days) ; 

1/rjh mean viremic period (in days); 

1/Hm average lifespan of adult mosquitoes (in days); 

if number of eggs at each deposit per capita (per day); 

1/HA natural mortality of larvae (per day) ; 

rjA maturation rate from larvae to adult (per day); 

m female mosquitoes per human; 

k number of larvae per human. 

The Dengue epidemic is modeled by the following nonlinear time-varying state equa- 
tions: 

' d it(t)= n h N h - [B(3 mh ^ h + //fc) S h 
^(t) = Bp mh ^S h - ( Vh + n h )I h 

' = V (l - dbfc) ( S ™ + Im) ~ (TJA + HA + <u) A m ' ' ' 

^|p(i) = T]AA m — ^B PhmJ^ + Hm + c m^ S m 
k (lP (t) — R ' (3fi m j^S m (/^m ~i~ Cm) Im 

with the initial conditions 

S/i(0) = Sho, hi®) = ho, -Rft(O) = Rho, ^ 
An(0) = A m Q, S m (0) = S m o, I m (0) = I m Q. 

3 Basic Reproduction Number 

An important measure of transmissibility of the disease is given by the basic reproduction 
number. It represents the expected number of secondary cases produced in a completed sus- 
ceptible population, by a typical infected individual during its entire period of infectiousness 

IH- 
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Theorem 3.1. The basic reproduction number IZq associated to the differential system ([!]) 
is given by 



n 



akB 2 /3 hm p mh M 

<p(Vh + Hh){Cm + Urn) 2 



Proof. The proof of this theorem is given in [9] . □ 

The model has two different populations (host and vector) and the expected basic 
reproduction number should reflect the infection transmitted from host to vector and vice- 
versa. If IZq < 1, then, on average, an infected individual produces less than one new 
infected individual over the course of its infectious period, and the disease cannot grow. 
Conversely, if TZq > 1, then each individual infects more than one person, and the disease 
can invade the population. 

Assuming that parameters are fixed, the threshold TZq is influenceable by the control 
values. Figure [T] gives this relationship. It is possible to realize that the control c m is the 
one that most influences the basic reproduction number to stay below unit. Besides, the 
control in the aquatic phase alone is not enough to maintain TZq below unit: it requires an 
application close to 100%. 



4 Optimal Control Approach 

Epidemiological models may give some basic guidelines for public health practitioners, com- 
paring the effectiveness of different potential management strategies. 
A cost functional was defined, 

ft/ r , 

(t) 2 + lL c A {t) 2 + IE (1 - a) 2 \dt, (3) 
Jo L J 

where jd, Is, 1l and 7^ are weights related to the costs of the disease, adulticide, larvicide 
and mechanical control, respectively. In this way, an OC problem is defined: 

minimize ( 3 ) 
subject to (1) and u2\. 



5 Numerical Implementation 

The simulations were carried out using the following numerical values: A^,, = 480000, B = 
0.8, p mh = 0.375, /3 hm = 0.375, fi h = 1/(71 x 365), i] h = 1/3, /i m = 1/10 )¥ > = 6, HA = 1/4, 
T]a = 0.08, m = 3, k = 3. The initial conditions for the problem were: S^o = — 10, 
I h o = 10, R h0 = 0, A m0 = kN h , S m0 = mN h , I m0 = 0. 
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(a) IZo as a function of c m and ca 




(b) IZo as a function of c m and a 




(c) IZo as a function of ca and a 
Figure 1: Influence of the controls on the basic reproduction number TZq 
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In a first approach, the same weights were considered, which means 70 = 75 = 7l = 
7£ = 0.25. 

The OC problem was solved using two different softwares: DOTcvp and Muscod-II. The 
simulation behavior is similar, and we decide to show here only the results of the DOTcvp. 

The optimal functions for the controls are given in Figure [2] 
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Figure 2: Optimal control functions (773 = 75 = = IE = 0.25) 

As the results for the basic reproduction number, the adulticide was the control that 
more influences the decreasing of that ratio, and as consequence the decreasing of the 
number of infected persons and mosquitoes. Therefore, the adulticide was almost the one 
to be used. We believe that the other controls do not assume an important role in the 
epidemic episode, because all the events happen in a short period of time, which means 
that adulticide has more impact. However the control of the mosquito in the aquatic phase 
can not be neglected. In situations of longer epidemic episodes or even in an endemic 
situation, the larval control represents an important tool. 

Figure [3] presents the number of infected human. Comparing the optimal control case 
with a situation with no control, the number of infected people decreased considerably. 
Besides, in the situation where optimal control is used, the peak of infected people is minor, 
which facilitates the work in Health centers, because they can provide a better medical 
monitoring. 
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Figure 3: Comparison of infected individuals under an optimal control situation and no 
controls. 



6 Conclusions 



A compartmental epidemiological model for Dengue disease was presented, composed by a 
set of differential equations. Simulations based on clean-up campaigns to remove the vector 
breeding sites, and also simulations on the application of insecticides (larvicide and adulti- 
cide), were made. It was shown that even with a low, although continuous, index of control 
over the time, the results are surprisingly positive. The adulticide was the most effective 
control, from the fact that with a low percentage of insecticide, the basic reproduction 
number is kept below unit and the infected humans was smaller. 

However, to bet only in adulticide is a risky decision. In some countries, such as Mexico 
and Brazil, the prolonged use of adulticides has been increasing the mosquito tolerance 
capacity to the product or even they become completely resistant. In countries where 
Dengue is a permanent threat, governments must act with differentiated tools. It will be 
interesting to analyze these controls in an endemic region and with several outbreaks. We 
claim that the results will be quite different. 
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